Multiband quasiparticle interference in the topological insulator Cua,Bi2Te3. 
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We present angle resolved photoemission experiments and scanning tunneling spectroscopy results 
on the doped topological insulator Cuo.2Bi2Te3. Quasi particle interference (QPI) measurements, 
based on high resolution conductance maps of the local density of states show that there are three 
distinct energy windows for quasi particle scattering. Using a model Hamiltonian for this system 
two new scattering channels are identified: the first between the surface states and the conduction 
band and the second between conduction band states. The comparison of the QPI data with self- 
energies extracted from photoemission indicates that spin-orbit terms are relevant for surface state 
-> conduction band scattering. 



The discovery of topological insulators (TIs)[T]-[3] has 
opened a new route to spintronic devices, whereby the 
surface states taking part in the charge transport are 
undisturbed by non-magnetic impurity scattering. Soon 
after the discovery of the topological phase in BiSb^^ [4 , a 
new family of TIs was predicted and discovered: Bi2Se3, 
Bi2Te3 and Sb2Te3 [5l|6l. Recently, it has become ap- 
parent that higher order terms in the k • p expansion are 
relevant for describing the low energy electronic states, 
leading to a warping of the Dirac cone away from the 
Dirac point [7, 8 . Such warping has been observed in 
angle-resolved photoemission spectroscopy (ARPES) ex- 
periments [9 and has been shown to be a source for local 
density of states (LDOS) oscillations observed in scan- 
ning tunneling microscopy /spectroscopy (STM/S) In 
terms of device applications it is important to understand 
and characterize such potential scattering channels, in 
particular in systems in which chemical potential tuning 
has taken place via bulk doping [9] . 

In this Letter we use STS, supplemented by ARPES 
and theoretical modeling, to uncover three different 
regimes of carrier scattering in Cu doped Bi2Te3. Firstly, 
scattering between the surface state (SS) and the bulk 
conduction band (CB) is shown to be vital in the un- 
derstanding of such doped systems. Secondly, warping- 
assisted SS-SS scattering is also observed and thirdly 
quasi particle scattering between CB states completes the 
line-up. 

Single crystals of Cuo.2Bi2Te3 were grown using a 
Bridgman technique. The elements were sealed in a 
quartz tube, heated to 850 °C and held there for 10 hours. 
This was followed by a slow cool (1 °C/hour) to 600 °C, 
the samples being held at this temperature for a further 
60 hours. ARPES was performed at the 1^ endstation 
on the UE112-PGM-2b beamline at BESSYII and with a 
lab-based system equipped with a high intensity helium 
discharge source. Angular and energy resolution were set 
to 0.002(0.008) A"^ and 4(15) meV at 1^ (lab-system). 
ARPES data were taken at 1^ at a temperature of 800 



mK, and the lab-based data at 18 K. UHV STM mea- 
surements were performed at 4.2 K, using Pt/Ir tips to 
achieve a tunneling current of 0.5 nA with junction resis- 
tances between 0.1-0.5 GVt. Conductance maps, ^(r, ^), 
have been recorded over 60 x 60 nm^ areas using a mod- 
ulation frequency of 969 Hz with a bias modulation am- 
plitude of 5 meV. In all experiments samples are cleaved 
in situ at room temperature. 

A powerful probe of quasiparticle scattering on 
TI surfaces is provided by conductance mapping in 




FIG. 1: (Color online) (a,b): LDOS maps measured at Vhias 
indicated. Insets show FT-LDOS, which have been three- 
fold symmetrized, (c): Dispersion of QPI ^-vectors along the 
V ^ K direction obtained from the FT-LDOS. Three different 
scattering regimes are indicated: I, SS^SS scattering due 
to Dirac cone warping; II, SS^CB and III: CB intraband 
scattering. 
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spectroscopic-imaging STM (TDHIS]. In Fig. [TJa,b) we 

present representative conductance maps at two differ- 
ent energies. There are clear standing wave-like patterns 
in g{r^E = 100 meV) originating at impurity sites, whose 
characteristic wavelength 30 A) can be picked up in 
the Fourier Transform (FT), p(q, cj), shown in the in- 
set. p(q, cj) at = -175 meV pTO] also shows a simi- 
lar FT pattern, whereby the lobes arranged with six-fold 
symmetry are clearly closer to the origin, making the 
associated longer wavelength patterns difficult to spot 
with the naked eye in the raw LDOS data. Many such 
LDOS maps have been measured, generating p(q, u) im- 
ages spanning the bias interval -350 meV to 250 meV in 
12.5 meV steps. At the energies where peaks are observed 
in p(q, cj), they always appear along the same directions 
in g'-space (F K). We determine the g'-vectors for the 
six quasiparticle interference (QPI) peaks and plot the 
average g'- vector in Fig. [T]^. Error bars are determined 
from the spread in the peak positions in the six direc- 
tions. We can identify three different regimes in the QPI 
patterns: (I) -300 meV to -175 meV, (II) -125 meV to 
125 meV and (III) > 200 meV. To understand the ori- 
gin of these different regimes, we compare the measured 
dl/dV to ARPES measurements presented in Fig. [2|a- 
d). Figure [2^ shows the Fermi surface (FS) of as-grown 
Cuo.2Bi2Te3. We observe a small hexagonal deviation 
(warping) from the perfect Dirac cone [71 [8] . A cut along 
r ^ M is shown in Fig. [SJd. Compared to pure Bi2Te3, 
Cu doping is expected to be p-type and we indeed find 
that the chemical potential is shifted below the bottom 
of the CB. From Fig. the energy of the Dirac point is 
found to be, Ejj = -180 meV and the valence band max- 
imum, Eyb = -80 meV. Exposure to ambient conditions 
is known to lead to n-type doping of Bi2Se3 [III [18], with- 
out destroying the topologically protected surface states. 
Figure [2]3,d shows that the same holds true for Cu-doped 
Bi2Te3 stored under ambient conditions prior to cleav- 
age: not only is the FS area clearly greater than for the 
pristine crystals, but the bulk CB states are also clearly 
visible inside the SS features, both in the FS-map and the 
I(k,£^)-image. Naturally, this n-type 'ambient doping' 
shifts the Dirac-point to higher binding energy. Fitting 
the T ^ K dispersion using a simple linear relation shows 
the Dirac point in the ambient-doped system to lie at -385 
meV, with a band velocity Vk ~ 2.52 eVA. This veloc- 
ity matches that in Fig. (2)3, further underpinning that 
the ambient exposure of the crystal does not significantly 
alter the properties of the SS. 

Since the STM conductance maps were also recorded 
on a crystal cleaved after exposure to air, we can now 
relate the different regions observed in p(q, cj) in Fig. 
[TJc) to features of the bandstructure. Region I spans 
the energy window between the top of the valence band 
(VB) and the CB bottom. In this special energy interval 
only quasiparticle scattering between surface states (SS- 
SS) is possible, facilitated where relevant by the Dirac 
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FIG. 2: (Color online) (a): Fermi surface contour of as-grown 
Cuo.2Bi2Te3 (T=800 mK, hu = 58eV). (b): l{k,E) along 
r ^ M . (c,d) same as (a,b) but for an ambient doped 
crystal, measured using hu = 21.2eV. 

cone warping [7 . For energies close to the bottom of the 
conduction band the QPI patterns are suppressed and 
then a new QPI regime is encountered (dubbed region 
II), matching an energy window in which both SS and 
CB states are available for scattering. In the following, 
we will show that region II is indeed dominated by scat- 
tering between surface states and the bulk conduction 
band (SS-CB). Finally, the dispersion observed at ener- 
gies above the second jump in the g'- vector around -hi 50 
meV - indicated in Fig. [TJc) as region III - is due to 
scattering within the conduction bands. 

To underpin the assignment of the different quasi parti- 
cle interference regions, and in particular the new insight 
regarding the multi-band scattering underlying region II, 
we utilize an extension of an electronic structure model 
based upon the work of Ref. [8 and applied to simulate 
QPI in Refs. [I9l[20]. The electronic structure of Bi2Te3 
and Bi2Se3 can be described by the Hamiltonian 

+ i^ir3(^^ - ^k^kl) + R2T^{'ikykl - kl) (1) 

where I4 is the identity and Fi-s are Dirac matrices sat- 
isfying the Clifford algebra {Ti^Vj} = 25ij and the func- 
tions £(k) = Co + CiA:^+C2A:^, M = Mo + MiA:^ + M2fc^, 
B = Bq -\- B2kl and A = Aq -\- A2k'^^. For systems such as 
Bi2Se3 and Bi2Te3, the last two terms in Eq. [l] break the 
in-plane symmetry, producing the Dirac cone warping. 



3 




D-cone D-cone+CB 

FIG. 3: (Color online)(a): Comparison between peccp(q, cj) and pcaic{^.,(^) from a simple Dirac cone model (left) and a Dirac 
cone model with the CB included (right) for energies of -175 meV and 50 meV. In between the pcaic{ci,(^) data are shown the 
relevant constant energy contours of the spectral functions with arrows indicating the dominant scattering processes, (b, c): 
QPI patterns for the D-cone+CB model at selected energies: -137.5 meV, falling in between regions I and II in Fig. [ijc), and 
at 213 meV where intraband scattering in the CB dominates (region III). 



This Hamiltonian can be reduced to an effective sur- 
face Hamiltonian, H2D = E{k) -\- Vk{kxCry — kyax) -\- 
|(/c^ + k^)crzj where and A are determined by com- 
parison to ARPES experiments. The latter Hamilto- 
nian has been used to study QPI patterns in STS ex- 
periments ^3 ^ and is able to describe SS^SS scat- 
tering. To include scattering between SS and the CB 
we need to calculate the surface Green's function start- 
ing from the model Hamiltonian, Eq. ([T]), using the 
method outlined in [21 [21]. The Hamiltonian is split 
into an in-plane component depending only on k\\ and 
a component depending on the interlayer hopping pa- 
rameters. The surface Green's function is calculated it- 
eratively whereby the Green's function for N layers, C^^, 
is related to the Green's function of a single layer, Go, 
through G^^ = Gq ^ — HIGn-iHz. Employing the fast 
iterative scheme proposed in [22 , we use the resulting 
surface Green's function to calculate p{q,uj) by means 
of the standard T-matrix formalism assuming ordinary 
impurity scattering with a momentum independent scat- 
tering potential [191 [20], 

To highlight differences between the QPI patterns ob- 
tained from the full, multiband Hamiltonian (labelled 
'D-cone+CB'), Eq. ([T]), and the effective surface Hamil- 
tonian (labelled 'D-cone'), we show results obtained for 
both models in Fig. [3^. In all calculations we used 



parameters chosen to reproduce the salient features of 
our ARPES data [23], and we added a chemical poten- 
tial to fix the position of the Dirac point to match the 
ambient-doped ARPES data. These ingredients allow 
the description of the multiband QPI and ARPES data 
in a self-consistent manner, and Fig. [3^ shows a com- 
parison of Pca/c(q, ^) and the STS results. For a bias 
voltage of -175 meV (top panel), both simulations show 
six dominant scattering vectors along T — K from intra 
SS scattering as indicated by the double-headed arrow 
in the constant energy contour, agreeing with the exper- 
imentally observed QPI in region I. The bottom panel 
shows the experimental data at a bias of 50 meV (re- 
gion II) compared to both models. The simple D-cone 
model predicts dominant scattering along the F — M di- 
rection, placing the QPI hexagon on its point, contrary to 
what is observed in the experimental p(q, cj). The main 
theoretical result is that this discrepancy is removed in 
the D-cone+CB model in which the outermost, intense 
features now match the orientation of the experimental 
QPI intensity and lie along F — as highlighted by the 
white arrow. The D-cone+CB model is also validated 
by the experimental QPI patterns at other energies. In 
Fig. the comparison is made for an energy in between 
regions I and II (top panel, -137.5 meV) and panel (c) 
for energy in region HI (213 meV). In the former, the 
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FIG. 4: (Color online) (a) ARPES dispersions (left) along 
r (red dots) and F (blue squares). The inset shows 
the expectations from k • p-theory. lmE{uj) from ARPES for 
the two directions is shown in the right-hand panel, (b) Angu- 
lar dependence of lmE{(jj) for two energies: purple circles for 
non- warped SS energy contours and black squares for warped. 
The energies at which the two lm^{uj) datasets are taken are 
indicated with dashed lines in (a). 

QPI is correctly predicted to be dominated by very small 
g'-scattering due to the very large DOS associated with 
the bottom of the CB. In region III, the Dirac cone+CB 
model suggests that the dominant scattering processes 
seen in experiment do not involve the SS at all, and are 
of CB^CB character. 

Fig. [3] shows we are able to successfully capture the 
essence of the experimental QPI behaviour by recogniz- 
ing that there are three scattering regimes. Based on 
this insight, we close the loop by extracting the g'- vectors 
for the dominant scattering processes from the two mod- 
els and overlay them as solid lines on the STS results in 
Fig. It is clear that the D-cone+CB model correctly 
explains both the size of the g'- vectors and the jumps be- 
tween the scattering regimes, whereas the single D-cone 
model is only able to describe region I. Closer inspection 
of Fig. ^ shows that the predicted energy dependence of 
the g'-vectors in the D-cone-hCB model does not match 
the experimental data perfectly. This points towards lim- 
itations of the k ' j9-theory in providing a sufficiently ac- 
curate description of the quasi-particle dispersions away 
from the Dirac point. Figure |4^ shows the experimental 
dispersion obtained from Lorentzian fits to the ARPES 
MDC's. The dispersion obtained within k • p-theory is 
shown in the inset of the same panel. While the warping 
in the experimental data appears as a reduced band ve- 
locity along F — M, k-p-theoiy predicts that the warping 
leads to an increase in velocity along T — K. Density 
functional theory is able to provide a more accurate de- 
scription of the true situation, correctly reproducing a 
lowering of the velocity along F — M upon the onset of 
warping [24 . 

A final point we address is the importance of spin-orbit 
scattering. The effect of spin-orbit terms in QPI scatter- 
ing has been studied in ref. [20 . Overall, spin-orbit 
scattering is not needed to capture the salient features 
of the QPI data. However, the angle dependence of the 
ARPES line- widths suggests that spin-orbit scattering is 



relevant in these systems. The right hand panel in Fig. 
|4^ shows the energy dependence of the imaginary part 
of the self energy, Im(i;(co')), obtained from ARPES and 
the local band- velocity. We note that the SS self energy 
is minimal (and close to azimuthal-angle independent) 
for energies above the VB maxima and below the on- 
set of warping, indicated by the purple points in Fig[4]3. 
Once warping kicks in, Im(S(co')) increases and becomes 
angle dependent (Fig. [4]3, black squares). Particularly 
this last effect calls for an angle-dependent interaction - 
the isotropic impurity scattering underpinning the QPI 
analysis will not do the trick - and the spin-orbit inter- 
action would provide just such an anisotropic scattering 
mechanism [20 . 
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